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Abstract. This paper describes an algorithm which produces all ideals compatible 
with a given surjective Frobenius near-splitting. 



1. Introduction 

Suppose that X is an algebraic variety over a field of characteristic p > 0. If X is 
Frobenius split with a Frobenius splitting cf) : F^Ox — >• Ox-, then X satisfies numer- 
ous remarkable properties. The following varieties possess a Frobenius splitting: toric 
varieties, Schubert varieties, ordinary Abelian varieties and, at least when reduced to 
characteristic p ^ 0, Fano varieties. In this context, it is very natural to study the 
compatibly (/)-split subvarieties (Z C X whose ideal sheaf Iz satisfies (j){F^Iz) ^ Iz)- 
These special subvarieties play a fundamental role whenever Frobenius split varieties are 
studied (see [BK05j ). Recent independent work in |KM09j . and independently by the 
second author, |Sch09j . have shown that there are only finitely many such subvarieties, 
also see [ShaOTj and |EH08j . In this paper, building on the ideas from |KM091[Sch09j . as 
well as ideas coming from tight closure theory, we exhibit an algorithm which computes 
all the compatibly (/>-split subvarieties. 

While Frobenius split varieties need not be affine, in this paper we restrict ourselves to 
affine varieties. This is not a terribly restrictive hypothesis since compatibly split subva- 
rieties of a projective variety X can be studied either on affine charts or by considering 
the affine cone over X. Our main result is as follows: 

Given a ring R and a surjective map (p : F^R — )■ R (for example a 
Frobenius splitting), we exhibit an algorithm which produces all the (p- 
compatible ideals. 

At each step, the algorithm produces the unique smallest non-zero (/)-compatible ideal, 
the so-called test ideal. 

Finally, we also explore a variant of this algorithm under the hypothesis that (j) is not 
necessarily a Frobenius splitting (or even surjective). This algorithm, and the original, 
have been implemented in Macaulay2 [KSllj . 
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2. Notation and background 

Convention 2.1. Through this paper all rings are commutative and of finite type over 
a perfect field k of characteristic p > 0, or they are a localization of such a ring. 

The algorithm in this paper indeed also holds on a much larger class of rings: at 
a basic level, this algorithm works for any ring of equal characteristic p such that the 
Frobenius map is a finite map (this condition is often called being F -finite). We use 
the notation (respectively F^R) to denote R viewed as an i?-module via Frobenius 
(respectively, via e- iterated Frobenius). More generally, for any i?-module M, F^M 
denotes the module M with the induced i2-module structure via Frobenius. Additionally, 
given an element r £ R, we will use F^r to denote the corresponding element of F^R. 
Finally, if / = (/i, . . . , fm) C is an ideal, we use /[P'l to denote the ideal {ff , . . . , /f ). 
This is easily seen to be independent of the choice of generators of /. 

Definition 2.2. We say that an /^-linear map cj) : F^R — t- i? is a splitting of (e-iterated) 
Frobenius, or simply a F -splitting, if (p sends F^l to 1. If i? has a Frobenius splitting, 
then we say that R is F -split (by (p). 

Definition 2.3. Given any i?-linear map <j) : F^R — ^ R (not necessarily a splitting), we 
say that an ideal J O R is (p- compatible if (f){F^J) C J, or simply that J is compatible 
with (j). If the (f) is clear, sometimes we will only say that J is compatible. If (j) is indeed 
a Frobenius splitting, then we say that J is compatibly ((p-Jsplit. 

Given (j) which is compatible with J as above, then there always exists a commutative 
diagram: 

(1) F^R > R 



f:{r/j)^r/j 

where the vertical arrows are the canonical surjections. We will use (j)/J to denote the 
induced map F^{R/J) — )• R/J as pictured above. 

The following well-known Lemma, which we will rely on heavily, follows immediately 
from the diagram above. 

Lemma 2.4. Assuming a commutative diagram ^^ as above, the (j)- compatible ideals 
containing J are in bijective correspondence with the (j)/ J- compatible ideals of R/J. 
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The following theorem motivates the main question of this paper. Its proof motivates 
the method of the algorithm 

Theorem 2.5. |KM09j . |Sch09| If cp : F^R R as above is surjective, then there are 
finitely many (f)- compatible ideals. 

Proof. Because the proof is motivating, we give a rough sketch of it here. Reduce to the 
case that i? is a domain. Fix a divisor Z)^ on the normal locus oi X = SpeciZ which 
corresponds to (j) (in the usual sense of Probenius splittings, see Proposition I2.6r b) and 
|BK05j ). Then by Proposition I2.6f b) below, all of the (^-compatible ideals have support 
contained within 

Supp(L») U Sing(X) 

where Sing(X) is the non-regular locus oi X. A roughly equivalent statement from the 
tight-closure perspective is as follows: all of the non-zero (^-compatible ideals contain the 
big test ideal ti){R,4>) (the unique smallest non-zero (/)-compatible ideal). 

Regardless, take Xi to be the union of the (/)-compatible subvarieties (vanishing loci of 
the (^-compatible ideals) . Either of the previous observations imply that the closure of Xi 
is a proper closed subset oi X = Spec R. Then repeat the process replacing X = Spec R 
by Xi and apply Noetherian induction. □ 

Therefore, in order to turn this theorem into an algorithm, one needs a way to identify 
the union of all subvarieties compatible with a given splitting (j). Equivalently and more 
algebraically, one needs to identify the smallest (/i-compatible ideal (that ideal is called 
a "test ideal"). 

We list some basic properties of Probenius splittings that we will need in what follows. 

Proposition 2.6. Suppose that R is as in Convention \2.1\ and that X = Speci?. We 

will assume that (j) G Hom/}(i<'^i?, i?) is a non-zero element. Then: 

(a) We have an isomorphism HomR{F^R, R) ^ F^Jj^~^'^ ^ F^T{X, Oxiil-p'')Kx)) 
of F^R-modules. In particular, if R is Gorenstein then Hom/j(F^^i2, i?) is a lo- 
cally free rank-one F^R-module. If}iomji{F^R,R) is free as an F^R-module, we 
label a generator of this F^R-module by . 

(b) By (a), cj) corresponds to an effective divisor D^j, linearly equivalent to {\—p^)Kx- 
Furthermore, for every (j)- compatible ideal J , V{J) X is contained within the 
set (SuppD^) U (SingX). 

(c) If (j) is surjective, then the set of cj)- compatible ideals is a finite set of radical ideals 
closed under sum and primary decomposition. 

(d) Let (/)* denote the composition 

Fl^R Fi^-^)^R ...^f:R^ R. 

Then if J is cp- compatible, it is also cp^ -compatible for any t > 0. If cj) is surjective, 
then every (jf -compatible ideal is also (j)- compatible. 

(e) // (j) is surjective, R is a domain, and a 7^ is any ideal that vanishes on the set 
(SuppZ),^) U (SingX), then 

4>{F:a)<Zcp\F^'a)<Z^^{F^-a)C... 

stabilizes to be the unique smallest (p- compatible ideal. 
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Proof. Parts (a) and the first part of (b) are in |BK05l Chapter 1] among many other 
places. The last part of (b) can be found in |KM09j but also is an immediately con- 
sequence of the theory of (sharp) test elements for pairs, see for example |Sch09j . The 
finiteness of part (c), as mentioned before was independently obtained in |KM09j and 
[SchODj . The other parts of (c) are the same as for compatibly split ideals and can be 
found in |BK05l Chapter 1]. 

The first statement in (d) is obvious. The second statement can be found in [SchlOl 
Proposition 4.1] although we also include a short proof for the convenience of the reader. 
To this end, suppose that 4f{Fl^J) C J for some surjective (j) ■ F*R — ^ R- It follows 
that is also surjective from which it follows that J is radical. All of the associated 
primes of J are also (/>*-compatible and since an intersection of (/)-compatible ideals is 
again compatible, it is harmless to assume that J is prime. Finally, a prime ideal is 0* 
or (/>-compatible if and only if it remains compatible after localization at itself, and so we 
may assume that J = m is a maximal (/>*-compatible ideal in a local ring (i?, m). But now 
suppose that (j){F^J) ^ J = m. In other words that the ideal 4){F^J) is not contained 
in the maximal ideal of a local ring, (i?,m). Thus (j){F^J) = R which certainly implies 
that (j)*{F^^J) = R and completes the proof of (c). 

We also prove (e) in some detail. First notice that o is necessarily contained in every 
non-zero ^-compatible ideal (because they are radical). Since (j) is surjective, there exists 
c £ R such that (j){F^c) = 1. For any z £ a, cz^ S a as well and so z = (t){F^zP c) G 
4){F^a). Thus a C (j)[F^a). Repeating this argument yields the ascending chain in (e). By 
the Noetherian hypothesis, this stabilizes say at C . By construction, it is a (/)-compatible 
ideal. If J 7^ is any other ^-compatible ideal then J = \[j 5 a by (b) and so the fact 
that C is minimal follows immediately. □ 

We point out a particular special case of (b) that we will use later. 

Corollary 2.7. Using the notation from Proposition \2. 6\( b). ifiiomR{F^R,R) is a free 
F^R-module with generator then (f)[F^») = ^ii{F^z-») for some z £ R. In this case 
= div(z). 

Consider the following example with regards to Proposition I2.6r a) . 

Example 2.8. Set S = k[xi, . . . , Xn] where A; is a perfect field of characteristic p. In this 
case, F^S is a free S-module with basis {F^x^^ . . . x^"}o<Xi<p^'-i- The map from 
Proposition I2.6r a) is the map which sends the basis element F^x^ ~^ . . . Xn ~^ to 1 and 
all the other basis elements to zero. We will explain later in Section [6] how such maps 
can easily be implemented in a computer. 

In this case, the divisor associated to $5- via Proposition I2.6r b) is the trivial divisor. 
Thus there are no non-trivial <I>5-compatible ideals by Proposition 12.6( b). 

We our next goal is to recall Fedder's Lemma; this will allow us to translate the 
problem of finding compatible ideals of R = k[xi, . . . , Xn]/I to finding compatible ideals 
on = k[xi, . . . , Xn]- The point is that if ii = S/I, then maps (f) : F^R — )• R come from 
maps (/) ■ F^S — )• S, which Fedder's Lemma precisely identifies. On the other hand, since 
$5 generates Hom5(F^^5, S), we can write ^s{F*z ■ •) = (f){F^»). Thus the choice of (j) 
is determined by the choice of a certain z £ S. Roughly speaking, Fedder's Lemma says 
that the set of allowable z are exactly /[f'l : /. 
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Lemma 2.9 (Fedder's Lemma). |Fed831 Lemma 1.6] Suppose that S = k[xi, . . . ,x„] for 
some perfect field k and R = S/I for some ideal I ^ S where p : S ^ R is the canonical 
surjection. Then: 

(a) If(j) : F^R R is any R-linear map, then there exists a S -linear map (j) : F^S — >• 
S which is compatible with I such that (p = (making the diagram commute 
as in EquationUl). 

(b) Given (j) and (p as in (a), then (j) = (p/I is surjective if and only if (j) is surjective 
at all points in a neighborhood of V{I) C Specs'. Furthermore, if 4>{F^») = 
^s{F^z ■ •), then (p is surjective at a point m G Spec 5" if and only if z ^ m'^"]. 

(c) An arbitrary map (p S Hom5(F^^5, 5) satisfies (p{F^I) I if and only if there 
exists c G /[P'] : / such that we can write (p{F^») = ^s{F^c ■ •) where $5 is as 
in Example \2.8[ Combining this with (a), we see that 

{F^{I^P'^:I))-Roms{F:S,S) 

is exactly the set of elements of Horns {F^S, S) which are compatible with /. 

(d) With the notation from (c), there exists an isomorphism: 

Yio^RiF^R, R) - : /)) • Homs(F,^5, 5)) / [{Ftl^^'^) ■ Hoxns{FtS, S)) . 

Proof. The proof found in jFed831 Lemma 1.6] is quite easy to read and so we will not 
repeat it here. □ 

We now state an easy corollary of Fedder's Lemma which shows that in fact we may 
reduce to the case of a Frobenius splitting on S — . . . , s^n]- 

Corollary 2.10. Using the notation from Lemma \2.W a). additionally suppose that <p is 
surjective. Then there exists a map ip : F^S — )• S such that 

(i) ip is a Frobenius splitting. 

(ii) For every cp-compatible ideal J R, the inverse image p~^{J) <^ S is compatibly 
split by tp (although there may be new tp-compatible ideals not coming from (p). 

Proof. For simplicity, we use M to denote the F^^S-module Hora.s{F^S, S). Fix any cp as 
in Lemma [2^91 a) and suppose that c G 5 is such that (p{F^c) ^ /, such a c exists by the 
surjectivity of (p. Now consider the -F^S-submodule W Q M generated by (F^I^p"^) ■ M 
and (p, in other words 

W:= {cP)Fes+{FtI^^'^)-M. 

Note that the formation of W commutes with localization. Also note that for any 
= (p{F^b ■ •)_+ ^s{FU ■•)eW, the induced map ip/I =: ijj on R = S/I simply 
coincides with (p{F^h ■ •) 

Consider now the map V :W ^ S defined by the rule V[ip) = ip{F^l). We will prove 
that r is surjective and thus that 1 G Image(r). It is enough to prove the statement 
locally at every prime Q £ Spec 5. There are two cases, either Q £ V{I) C SpecS or 
not. 

(a) If Q G then after localizing at Q, we notice that cpQ : F^Sq — )■ Sq has 1 
in its image, say (pQ{a) = 1. It follows that the function 7(F^^») = (pQ{F^a ■ •) 
sends 1 to 1 and so rQ(7) = 1. 

(b) If Q ^ ^(-^)) then it is even easier since then Wq = Hom5Q(F^S'Q, Sq) and claim 
follows immediately. 
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Since now 1 G Image(r), there exists ip S }ioms{F^S, S) such that ip{F^l) = 1. 

Suppose finally that J C i? is i^-compatible. We must show that J' = p~^{J) is tp- 
compatible. We notice that ^p induces tp/I : F^R — )• R and furthermore, tp/I{F^») = 
4>{F^b ■ •) for some 6 G i? as observed initially. It is enough to show that tp/I{F^J) C J 
but this is obvious by the above characterization oi tp/I. □ 

As a final remark in this section, we mention a common source of compatibly Frobenius 
split ideals. 

Remark 2.11. Suppose that X is a projective Frobenius split variety projectively nor- 
mally embedded into P". In that case the affine cone over X is Frobenius split and 
so yields a Frobenius split ring (the projective normality then guarantees that it is a 
quotient of the ring S = k[xQ, . . . , Xn\)- This is a particularly common way of producing 
Frobenius split rings. In this case, the same Fedder-type Lemma works, and one can 
compute the compatibly Frobenius split subvarieties of X by computing the compatible 
ideals of Spec 5. 

3. The statement of the algorithm 

Suppose that S := A;[xi, . . . , a;„,] where /c is a perfect field of characteristic p > 0. 

Because of Fedder's Lemma 12.91 and Corollary 12.101 we additionally suppose that 
(/> : F^S —7- 5 is a surjective 5'-linear map (for example a Frobenius splitting) that is 
compatible with / (and henceforth / will not play much of a role). In a later section, we 
will handle the non-surjective case. The advantage with working with S instead of R is 
that F^S is a free 5-module, and so specifying (p is the same as specifying where a basis 
is sent. Finally we fix z G 5 such that ^s{F^z ■ •) = (j){F^») where $5 is as in Example 
[2:51 

With the notation above, and given any prime ideal Q ^ S {va practice containing 
/) the recursive algorithm described below produces a list of all prime i;^-compatible 
ideals which properly contain Q. In order to do this, the algorithm finds the smallest 
0-compatible ideal properly containing Q. The initial input Q to the algorithm can be 
the zero ideal or each of the minimal primes of /. 

Of course, the plan of the algorithm is to apply Proposition I2.6r e) to the ring S/Q, 
see Step (3) below. Therefore we need to define the ideal a which all the Frobenius split 
subvarieties contain. Identifying the singular locus is straightforward, see step (1) below, 
and identifying the locus corresponding to the divisor on R is accomplished in (2). 
Step (4) is then recursive where we replace Q by larger ideals. 

Here are the steps of the algorithm (later we will describe an algorithm which works 
for non-surjective (p). 

(1) Find an ideal J S such that Spec(5/Q) \ {V{J) R Spec(S'/Q)) is a regular 
non-empty scheme. For example J could define the singular locus of S/Q. 

(2) Compute B := Ann^ {{Q^'^ : Q)/{{z) + QIp^I)) = {{{z) + Q^p^^) : (QIp^I : Q)) . 
This ideal in not contained in Q as long as (p/Q is non-zero. 

(3) Find the first t such that (p^F^^JB + Q)) = cP^+^f!;^^'^^^ {J B + Q)) = C (as in 
Proposition 12.6( e)). We will show that the ideal C properly contains Q. 

(4) Let Qi, . . . , Qm be the minimal primes of C, add them to the list of compatible 
ideals, and repeat the algorithm with Qi = Q. Every i;^>-compatible ideal, properly 
containing Q, that is minimal with respect to inclusion, appears in the list of the 

Qi. 
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In Section [6] below, we will discuss issues of complexity, especially with regards to the 
integer t needed in Step (3), 

4. The proof that the algorithm works 

In this section, we prove the individual claims made in the algorithm. Step (1) has 
no associated claims so we move on to Step (2) which is the only technical point in the 
argument. Let us briefly explain the difficulty. While the divisor associated to (j) is easy 
to identify on SpecS* (it is merely div(z)), it is harder to identify on S/Q (which need 
not even be normal), indeed, this divisor has no transparent relation to div(2;) except in 
special cases. Step (2) identifies an ideal whose support contains -D^/q at least on the 
locus where S/Q is normal. To do this we utilize Fedder's Lemma 12.91 

Let us informally explain where the ideal B comes from. Indeed, (Qt^'l : Q) cor- 
responds via Fedder's Lemma to the set of all elements of Iioms{F^S, S) which are 
compatible with Q. On the other hand, (z) + Q^^"^ corresponds to the submodule of 
HomsiF^S, S) generated by all tp which are both compatible with Q and such that 
ip/Q = (p/Q- We then observe that 

(((z)+Q[*'^l):(Q[f^]:Q)) 

simply defines the locus where these two modules are distinct. On the regular locus of 
R/Q, this is simply the defining equation of D^/q. 

Now we carefully prove a generalization of the claim from step (2). 

Lemma 4.1. For any (p '■ F^S — S compatible with Q, the ideal B constructed in (2) is 
not contained in Q as long as (j){F^S) is not contained in Q (we do not require that (p is 
surjective). 

Proof. Consider the induced S/Q-linear map (p/Q : F^{S/Q) — )■ S/Q. Notice that 
Image((/>/Q) = (p{F^S)/Q / since Q does not contain (p{F^S). Therefore, the map 
cp/Q : F^{S/Q) — )■ S/Q is not the zero map. We will use N to denote the F^^5-module 
Hom5(F,"5, S) and use M to denote the F,^(5/Q)-module YIouis/q{F^{S /Q), S/Q). Con- 
sider the cyclic F^^(5/Q)-submodule of M 

K := {(P/Q) C M := Hom5/Q(F,^(5/Q), S/Q). 

Now, M is a F^[S/Q)-m.odvle of generic rank 1, and so there exists an element F^d G 
F^S \ F^Q such that 

{Ftd) -MCK. 

We will show that d G F^B which will imply that B is not contained in Q. 
We use Fedder's Lemma 12.91 to see that 

Boms/QiF:{S/Q), S/Q) = M - (f,^(Q[p^1 : Q) ■ n) / [{F^Q^p'^) ■ n) 
Furthermore, under this identification, the submodule K corresponds to 

[f^Hz) + QlP^l) • n) I ((F^QlP^l) • n) . 

It follows that d multiphes Fl{(Q^v'^ : Q) into F^^((z) -FQ[p'1) and so d £ B as desired. □ 
The remaining claims from the algorithm are easy and we now prove the claim in (3). 
Lemma 4.2. The ideal C constructed in (3) exists and properly contains Q. 
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Proof. By (1) and (2), the ideal JB is not contained in Q, thus JB + Q properly contains 
Q. The result then follows since this chain of ideals stabilizes by Proposition I2.6r e) . □ 

Finally, we prove the claim in step (4). 

Lemma 4.3. With the notation of step (4), if P is any prime (p- compatible ideal, properly 
containing Q, then P contains Qi for some 1 < i < m. 

Proof. This follows immediately from Proposition 12.6( e). □ 

5. A MORE GENERAL ALGORITHM 

In this section we describe a slight modification of this algorithm which produces 
something of interest even if (j) is not surjective. Fix S and <j) : F^S ^ S an S- linear 
map as in section [Sj but do not assume that (j) is surjective. Consider the ideal of 
S, K = Y^Image((/)). With a minor change to the algorithm presented above, we can 
compute all the prime i;^-compatible ideals of J not containing K. Furthermore, we can 
identify all the (/^-compatible primes, see Remark 15.31 

Obviously there is another way to do this too, set V{K) C SpecS to be the vanishing 
locus of K. Then one can cover Speci? \ V{K) with affine charts, charts where (j) is 
surjective, and compute the compatible ideals on each of those charts. However, a minor 
change of our original algorithm allows us to run it without this complication. 

The input for each stage of our algorithm is the same as before and the steps are quite 
similar (listed below): 

(1*) Find an ideal J C 5 such that Spec(5/Q) \ {V{J) n Spec(5/Q)) is a regular 
non-empty scheme. For example J could define the singular locus of S/Q. 

(2*) Compute B := Anns {{Q^"^^ : Q)/{{z) + Q^^^)) = {{z) + Q^^'^) : (Qt^'l : Q)) . 
This ideal in not contained in Q as long as Q does not contain K. 

(3*) Define an ascending chain of ideals of S recursively as follows: Co = JB + Q and 
Ct = (t){F^Ct-i) + Ct-i. Find the first t such that Ct = Ct+i and set C = Q. 
The ideal C properly contains Q. 

(4*) Let Qi, . . . ,Qm be the minimal primes of C which do not contain K. Then 
repeat the algorithm with Qi = Q (all of the Qi will be (/>-compatible ideals so 
should be added to the list of valid outputs). Every (^-compatible ideal, properly 
containing Q but not that is minimal with respect to inclusion, appears in 
the list of the Qi. 

Running this recursively will produce all (/)-compatible primes not containing K. 

Now we prove that the algorithm is correct. Again property (1*) has no associated 
claims. Property (2*)'s proof is already contained in the proof of Lemma l4. 11 We now 
prove that the assertions in properties (3*) and (4*) hold. 

Lemma 5.1. The ideal C defined in (3*) properly contains Q. 

Proof. This chain also is ascending and is thus eventually constant. It also properly 
contains Q since Co does. □ 

Lemma 5.2. With the notation of step (4*), if P is any prime (j)- compatible ideal, 
properly containing Q, then P contains Qi for some 1 < i < m. 
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Proof. Note that C vanishes on the same locus as the test ideal t{S/Q, (j)/Q) by |Sch09] 
and Proposition 12. 6[ In particular, P contains y^T{S/Q, (j)/Q) 5 {JB + Q)/Q = Cq 
because the test ideal is the unique smallest (;^/Q-compatible ideal. Thus P = (j){F^P) + 
P D (j){F^Co) + Cq = Ci and so recursively, P contains C and thus also P contains ^/C 
and so it also contains a minimal prime of ^/C. □ 

Remark 5.3. Finally we explain how to find all the i^-compatible primes. We have just 
found all the i;^>-compatible ideals not containing K. On the other hand, suppose P ^ K 
is a prime ideal, we will show it is always i;^>-compatible. It is easy to see that P is (j)- 
compatible if and only if P/K is (p/K compatible. But P/K is clearly compatible since 
(f>/K IS zero. 

6. Further remarks 

In this section we briefly discuss issues related to the implementation of this algorithm 
and also discuss some connections with previous work. 

6.1. Notes on implementation and complexity. In this subsection, we briefly dis- 
cuss the issues surrounding the implementation of this algorithm. 

The computationally intensive steps in the algorithm involve computation of the sin- 
gular locus in step (1), the computation of colon ideals in step (2), the repeated applica- 
tion of (j) in (3) and finally primary decomposition in step (4). Indeed, all the steps (1), 
(2) and (4) are already implemented in many computer algebra systems (for example 
Macaulay2). Step (1) is simply the computation of, potentially many, minors (although 
any single minor that doesn't vanish on the given Q would suffice). The computation 
of the colons of ideals in step (2) reduces to the computations of ideal intersections as 
described in |GP08l Section 1.8.8]. Step (3) will be handled below and the primary 
decomposition in step (4), a list of references for the history of computing primary de- 
composition can be found on |CLO071 Page 206]. 

In order to implement the algorithm (or the generalized version from Section [5]), we 
need to answer the following question. 

Question 6.1. How does one compute the images of the S'-linear maps (p : F^S ^ S7 

We fix a = (F» • $5 e Roms{F^S,S). Given any ideal J C 5, ^s{F^J) is an 
ideal that has appeared previously in several contexts. Again, is as in example 12.81 
Notably, in |BMS08j and |Kat08j . also see [BSTZlOl Proposition 3.10], it was shown that 
^siF'^J) is the unique smallest ideal A (Z S with the property that J C A^^h In these 
contexts it was denoted by jI^/p^I and Ie{J) respectively. In the context described, where 
5 is a free module over , this ideal is highly computable as we next show (cf. |Kat08t 
Section 5] and [BMSOSl Proposition 2.5]). 

It is easy to see that for ideals Ji, J2 C S, we have /e(Ji -|- J2) = Ie{Ji) + leiJi)-, so 
the calculation of Ie{J) reduces to the case where J is generated by one element g G S. 
Let a denote n-tuple of non-negative integers (ai, . . . ,a„), let < a < denote the 
condition < ai, . . . , a„ < p'^ and write Now write 

bee, 

0<a<p'= 

where gi,^a S S. We claim that leiid)) is the ideal A generated by {gb,a \ b £ B,0 < a < p^} 
To see this note first that, clearly, g G A^^"^. If L C S" is such that g G L^^"^ then we can 
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find ai, . . . ,as € L and ri, . . . , S S such that 



E 



6G8, j=l 
0<a<p<= 



For all 1 < z < s we can now write 



0<a<p<= 



where r;, „ j G S and we obtain 



E 

tee, 

0<a<p'= 



P 



E 

beB, 

0<Q<p = 



Since these are direct sums, we may compare coefficients and deduce that for all 6 G 5 
and 0<a<p^, g^^^ = ^^^i rfa^^af hence gb,a = I]i=i and gb,a e L. 

This construction translates easily into an algorithm as follows. Extend the ring S to 
T = S[yi, . . . , Un], chose a term ordering in which the variables bigger than 

yi, . . . ,yn and reduce the g £ T with respect to the ideal generated by —yi, ■ ■ ■ ,Xn —yn 
to obtain g = X]o<o<p<= where G /^[i/i, . . . , For each < a < p*^ we can write 

ga as a sum of terms Aiy'^^^' + • • • + XmU^''"^^ with Ai, . . . , Am G k and for each 1 < i < m 
we can write A, = ^^eB ^Ib^ where Aj,b G /c. Now ^q, = T.b(^B YaU K^b^v'^^'^ and /e(5') is 
the ideal generated by X^^^L^ Aj^fcx'^**^ for all choices of < a < p*^ and h £ B. Notice that 
the complexity of applying (f) is essentially the complexity of finding the S'^'-coordinates 
of elements in S in terms of a given set of free generators of S. 
One can also ask the following. 

Question 6.2. How many times must (f) be applied in step (3)? 

In the experiments we have done so far, the condition in step (3) does not seem to be 
a limiting factor since applying is itself a very fast operation and the required value 
of t is quite small. Of course, the particular t needed also depends upon the given ideal 
Q. 

Indeed, one can give a reasonable bound on t based upon the sort of analysis found in 
|BMS08j . Fix a an ideal in 5 = k[xi, . . . , x„] and write (t>{F^») = ^s{F^z ■ •). Suppose 
that a is generated by elements of degree at most d in 5 so that z • a is generated by 
elements of degree at most {deg z)+d. We now consider the generators of F^{z-a.) C F^S 
as an S -module. It is easy to see that F^{z ■ a) is generated by elements of degree at 
most (degz) + d + (p^ — l)n (since F^S is as an 5-module is generated by elements of 
degree at most (p*^ — l)n). Note that if / is a polynomial of degree k, then it follows that 
^s{F^f) is a polynomial of degree at most ^^j^-^^J. We the see that ^s{F^z ■ a) is 
generated by elements of degree at most 



(deg z) + d + {p*^ — l)n — {p^ — l)n 




(deg z) + d 


p^ 




p^ 
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If we apply (j) again, we obtain an ideal generated by elements of degree at most 



(deg z) + 



(deg2:)+d 



< 



(1 +p^)(deg2;) +d 



,2e 



For i » 0, (j) {Fl^a) is generated by elements of degree at most 



(1 + + • • • + p(^"^^)(deg z) + d 



< 



l^pte _ l)((ieg z) 



p 



{t+l)e 





deg z 


+ 1 < 









+ 1. 



We now ask how many times need we to apply (j) before we reach this stable degree 
(which is a vector space of polynomials of bounded degree). But for this we merely need 
d/p^^ < 1, or in other words after at most [logpe((i)] applications of cj). 

Of course, we still may need to apply (p further when are working within this vector 
space. However, as soon as the containment (j)^{F^'^z ■ a) 5 </)*"^^(i^i*'''^'*'^z • o) is equality, 
that step in our algorithm terminates. In particular, we are working within this fixed 
vector space which has dimension at most the binomial coefficient: 



M := 



dee 2 



+ 1 + n 
n 



Recall that n is the number of variables in S. Therefore, we require at most M applica- 
tions of 4>. 

In summary, in step (3), we need only apply (j) at most: 



riog„e(d)i + 



deg Z 



+ 1 + n 
n 



times. 



6.2. Connections with previous work. We now also briefly explain some of the other 
ways the ideas in this algorithm have previously appeared. We use the notation from 
Section [6711 in particular (/)(•) = ^s{u ■ •). 



Definition 6.3. |Kat081 Definition 5.5], c/. |BMS08j As above, fix e i?. For any ideal 
J C S we define J* to be the smallest ideal A ^ S containing J with the property 

uA C A^p^y 

It follows that the sequence of ideals {Ci}i>o from step (3*) stabilizes at the value 
{JB + QY This construction (and the fact that it had already been computed by the 
first author) was part of the inspiration for this project. In view of this, step (3*) can 
be totally phrased in the language of J* 

The original motivation for studying J* " was to compute the test ideal. In particular, 
the ideal C we construct in (3*) restricts to the test ideal T{S/Q,(j)/Q) in many cases 
(with additional work, it can always be made to restrict to the test ideal, see for example 
|Kat08j ): however, it always restricts to the test ideal T{S/Q,(j)/Q) if cp/Q is surjective. 

Remark 6.4. One can also view the results of this paper from the point of view of 
Frobenius maps on the injective hull of residue fields. Let {S, m) be a complete local 
regular ring and E = Es{S/m) is the injective hull of its residue field. Let / : 5 — )■ 5 
be the Frobenius map f{s) = s^, and let S[@; /^] be the skew-polynomial ring with 
coefficients in S where the variable O satisfies Qs = f^{s)Q for all s £ S. 
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The 5-module E has a natural structure of S[T; f^]-modu\e which can be described 
by identifying E with a module of inverse polynomials k[x^ , . . . , x~] (cf. |BS981 Example 
12.4.1]) and extending additively the action TXx^^^ ■ ■ ■ Xn^"' = A^^x^^ . . . ^" for 
all A € A; and /3i, . . . , /3„ > 0. Any structure of S[@; /^]-module on E is given by = uT, 
where T is the natural action above, and, with this G, an 5-submodule Ann^; J <^ E 
is an S[@; /'^]-submodule if and only if uJ C jtp"! (cf. [KatOSl section 4]). Thus we see 
that the uT-compatible ideals oi R = S/I are the annihilators of 5[6; /'^J-submodules 
of E which contain J, i.e., /'^J-submodules of Ann^;/ = En^R/mR) and hence 
our algorithm produces these. These annihilators form the set of special ideals in the 
language of |Sha07j and |Kat08| . An analysis of our algorithm shows that it will produce 
all special primes P for which the restriction of Q to Ann^; P is not the zero map. 

7. The algorithm in action 

In this section we present some interesting calculations performed with a Macaulay2 
implementation of the algorithms presented in this paper. 

First we include an example where we step through the algorithm. This example is 
also interesting because the ideal defining the singular locus of i? = S/I is not always 
compatible with our choice of (p. 

Example 7.1. Consider the ring S = k[x, y, z, w] where k is any perfect field of charac- 
teristic 3 and set / = (x^ — yz). 
We set 

/ 2 \2 2/ ,i\ 62,4 2,2222,52,3 2, 222 

z = [X — yz) w x[x + 1) = X w + X yzw + x y z w + x w + x yzw + xy z w 

and fix 4>{F^,») = <ts{F^,z ■ •). It is easy to see that cj) is /-compatible by Fedder's 
Lemma l2 . 91 noting that 2; is a multiple of {x'^w — yzw)'^ . Furthermore, note that z has a 
term x'^y'^z'^w'^ ^ m'^] = {x^,y'^,z^,w'^), which implies that (f) is surjective at the origin. 
More generally, the same term implies that (/)(-F*l) = 1 which means that (/) is surjective 
everywhere as it is a Frobenius splitting. 

First we set Q = / In step (1) of the algorithm, we compute the singular locus of S/Q, 
it is defined by the ideal J = {x,y,z). For step (2) of the algorithm, Macaulay2 will 
easily verify that B := {z + Q^^^) ■ (Q'^' : Q) = {x'^ — yz, x'^w'^ + xw'^). Now we move on 
to step (3): one can either verify by hand, or by Macaulay2 that 

(j){F^JB + Q) = {w,x^ -yz). 

On the other hand 

02 {FliJB + Q)) = cl){F,cl){F,{w,x^ - yz))) = {w,x^ - yz) 

as well, so the ideal C = {w,x'^ — yz). This ideal is already prime so no primary 
decomposition is needed for step (4). 

Now we repeat the algorithm with a new Q' := {w,x'^ — yz). In step (1), the singular 
locus is now J' = {x, y, z, w). In step (2), one obtains B' := {w, yz + x,x'^ + x). For step 
(3), we easily compute that 

cPiF.J'B' + Q') = S = cl)\F^J'B' + Q') 

In particular, there are no new Qi and the algorithm terminates. 

It follows that the only proper non-zero i;^-compatible ideal of S properly containing 
/ is {w,x'^ — yz)s. In particular, the only proper non-zero cji/I = 0-compatible ideal of 
R = S/I is the ideal {w)ii. 
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Now we perform a more involved calculation which looks in greater detail at the 
example given in |Kat08[ section 9]. 

Example 7.2. Let K be the field of two elements, S = X2, X3, X4, X5], denote 

m = (xi,X2,X3,X4,X5). 

Now take / be the ideal of S generated by the 2x2 minors of 



X4 X4 



X2 X5 
X3 Xl 



and let S = R/I. We consider a 
i?-linear map from Example 

^ I ^3 , _2 



) G }iom.ji{F^ S , S) induced by pre- multiplying the 
by an element z G (jt^^ '■ I) S 

, _2_2_. , 2_2 , „3^2 



Z = X1X2X3 + X1X2X4, + X1X3X4X5 + X1X2X3X4X5 + XlX2X4^X5 + X2X4X5 + X3X4X5 + X^X^ 

2 
5 



Since z ^ (xf , Xg, 3;|, Xg) = m'^^, we see that (j) is surjective from Fedder's Lemma, 
Lemma 12.9( b). 

Our algorithm now produces a complete set of 0-compatible primes as follows: 

R, (xi,X4), (xi,X4,X5) 
{Xi + X2,xl + X4X5), {xi + X2,x\ + X4X5), (X3 + X4,Xi + X2.,x\ + X4X5), 
(X1,X2,X5,X3 +X4), (xi,X2,X4), (xi , X2 , X5) , (xi , X3, X4) , 
(xi,X2,X3,X4), (xi,X2,X4,X5), (xi , X3 , X4, X5) , m 

Consider now the Frobenius action = zT on the injective hull E of the residue 
field of ]K|xi, X2, X3, X4, X5]. The ^-compatible primes above are also the the special 
primes of the S[Q] /]-module E. In |Kat08( section 9] it was shown that there is a 
S[Q; /]-linear surjection of E onto H^g(S') where the latter is equipped with its canonical 
5[G; /]-module structure, and the set of prime annihilators of S[Q; /]-submodules of this 
quotient consist of the first three (/)-compatible primes above. 

Example 7.3. We now consider an example suggested by the referee and inspired by 
the calculation of Schubert varieties, see in particular |Knu091 section 7] for the origin 
of the element u below. 

Let S to be a polynomial ring in indeterminates {x^j | 1 < j < i < 4} over a field IC of 
prime characteristic 2 and let 

U = X4i(x3iX42 - X41X32)(X41 - X21X42 - X31X43 + X21X32X43); 

this is the product of the four lower left minors of the matrix 

10 

X21 1 

X31 X32 1 

X41 X42 X43 1 



M 



For sets a, /3 C {1, 2, 3, 4} of the same cardinality, let [a, /3] denote the determinant of 
the submatrix of M consisting of the rows in a and the columns in /3. Our algorithm 
produces the following twenty-three compatible ideals: the ideal generated by all vari- 
ables, twelve ideals generated by the variables in the positions 















i:i 


c: 
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iil M\ Sii 

together with the ten ideals 
([34, 12]), ([34, 12], [34, 13], [34, 23]), ([34, 12], [23, 12], [24, 12]), 
([34, 12], [23, 12], [34, 13], [34, 23], [24, 12]), ([234, 123]), ([234, 123], X41), 

([23, 12],X41,X42), ([23, 12],X41,X42,X43), 

( [34, 23] , X31 , X41) , ( [34, 23] , X21 , X31 ,Xii), 

Finally, we compute an example when the map (p is not surjective to illustrate the 
more general algorithm. 

Example 7.4. We fix a 2 x 4 matrix of indeterminates, and we let S* to be a polynomial 
ring in these indeterminates over a field K of prime characteristic 2. 

For any 1 < i < j < 4 we denote Ajj the 2x2 minor obtained from columns i and 
j and for any subset A C {1,2,3,4} we denote Va the ideal generated by the matrix 
entries in all columns listed in A. 

For z = A12A13A14 = (xiia;22 - X2ia;i2)(a;ii2;23 - a;2ia;i3)(a;iia;24 - X2iXii), we form 
</>(-F*») = ^^(-F*^-*). This (j) is easily seen to not be surjective. Our generalized algorithm 
produces two sets of compatible ideals. The first is the poset 



(Ai2,Ai3,Ai4,A23,A24,A34) 




It is easy to see that z and so the induced map 0/14 is the zero map since 

(f){F^,S) C In particular, our algorithm also produces the following poset of primes. 



^1,2,3,4 




Our generalized algorithm is only guaranteed to find those primes which are not contained 
in (j)^{FlS) for i ^ 0. However, these primes contain no information, see Remark 15.31 
and are simply an artifact of the algorithm. 
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